Preprint typeset in JHEP style - HYPER VERSION 

Holographic mesons in D4/D6 model revisited 



Kwanghyun Jo 

Department of Physics Hanyang University, Seoul 133-791, Korea 



E-mail: j okh38@gmail . com 



Youngman Kim 

Asia Pacific Center for Theoretical Physics and Department of Physics, Pohang 
University of Science and Technology, Pohang, Gyeongbuk 790-784, Korea 



E-mail: ykimOapctp . org 



Sang-Jin Sin 

Department of Physics Hanyang University, Seoul 133-791, Korea 



E-mail: 



sjsinOhanyang . ac .kr 



Abstract: We revisit holographic mesons in the D4/D6 model to study holographic light 
vector mesons and the properties of heavy quarkonium in confining and deconfining phases. 
To treat the light mesons and heavy quarkonium on the same footing, we use the same 
compactification scale Mkk in both systems. We observe that like scalar and pseudo-scalar 
mesons the vector meson mass is linearly proportional to the square root of the quark mass, 
when the quark mass is large. With a Mkk fixed by light meson masses, we calculate the 
mass of heavy quarkonium in confining and deconfining phases. We also obtain the in- 
medium dispersion relation for heavy quarkonium, which is important to understand the 
dissociation and the screening mass of heavy quarkonium in quark-gluon plasma. 
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1. Introduction 

The properties of the heavy quarkonium both at zero and at finite temperature have been 
intensively studied, see for a review. At zero temperature, the charmonium spectrum 
reflects detailed information about confinement and quark-antiquark potentials in QCD 0|. 
At finite temperature, melting of heavy quarkonia could be a signal of the formation of 
the quark-gluon plasma (QGP) in a relativistic heavy ion collision [||. Moreover, lattice 
calculations suggest that the charmonium states will survive at finite temperature up to 
about 1.6 to 2 times the critical temperature T c ||, |5|. This suggests that analyzing the 
charmonium data from heavy ion collision inevitably requires detailed information about 
the properties of charmonium states in QGP. Therefore, it is very important theoretical 
challenge to develop a consistent non-perturbative QCD picture for the heavy quark system 
both below and above the QCD phase transition temperature. In this respect, a promising 
attempt would be the holographic QCD (via AdS/CFT §§). 

In a stringy D4/D6 model, scalar and pseudo-scalar bound states have been studied 
at zero and finite temperature in (?], B. In a more phenomenological approach, bottom-up 
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model, the mass spectrum of a charmonium and its dissociation temperature have been 
investigated 10, plf| . However, we note that in the bottom-up model, different infrared 
scales are introduced to describe light mesons and heavy quarkonia. For example, in the 
hard wall model, the location of the infrared cutoff z m varies from light mesons to heavy 
quarkonia: l/z m ~ 320 MeV p|, 111 for the light meson and l/z m ~ 1315 MeV |9| for the 
charmonium. 

In this paper we study the spectrum of the light vector meson and heavy quarkonium 
using the D4/D6 model Q in high-temperature deconfining phase as well as in confining 
phase. To treat the light mesons and heavy quarkonium on the same footing, we use 
the same compactification scale Mkk in both systems. For the light vector meson, the 
spectrum was discussed in the D4/D8/-D8 model [14] where the chiral symmetry and its 
breaking are realized geometrically. However, in the Sakai-Sugimoto model [14] it is quite 
difficult to include the quark mass. Therefore it is still of worth to study the spectrum of 
the vector mesons in the D4/D6 model and study the effect of the quark mass. We observe 
that like scalar and pseudo-scalar mesons [j/J the vector meson mass is proportional to 
a square root of the quark mass M 2 ~ m q for large quark masses. We also obtain the 
in-medium dispersion relation for heavy quarkonium, which is important to understand 
the dissociation and the screening mass of heavy quarkonium Jl^| . 



2. D4/D6 model 

We briefly summarize a pioneering holographic QCD model, the D4/D6 system j7|. The 
model contains N c number of D4 branes and Nj flavor D6 branes whose configuration is 
given in Table 1. In the probe limit, the N c D4 branes are replaced by their supergravity 
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Table 1: The brane configuration : the background D4 and the probe D6 



background, and Nj D6 branes are treated as probes. In this model mesons of a QCD-like 
gauge theory are described by the fluctuations of the D6 brane in the D4 background. The 
geometry of confining D4 brane reads 

3/2 / j \ 3/2 



ds 2 = (j) ( V ^dx»dx» + f(U)dr 2 ) + (£\ (^L- + u^Q 



e^ = 5s (£) 3/ \ F, = dC 3 = ^e 4 , /(£/) = 1-^| , (2.1) 
where string coupling constant g s and the period of r are given by 

q 9 2 YM 6t =^Il!1 (22) 
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The parameter L is given by the string coupling constant g s and the string length l s , 
L 3 = Trg s N c lg, and the compactification scale Mkk reads 

By introducing K(p), the metric is simplified to be 

/TT\ 3/2 

<fc 2 = ( _ j V ^dx»dx u + K(p) [dA 2 + A 2 ^ 2 , + dr 2 + r 2 d</> 2 ] , (2.4) 

where 

7-7-1/2 / rr3 \ 2/3 

K(p)^L 3 / 2 ^, C /( /9 ) = ph+^J wherep 2 = A 2 + r 2 . (2.5) 

The position of the D6 brane is described by r(A) with (p=0 and r =constant. Then the 
induced metric on D6 is 

//T\3/2 

d4 6 = ( j J V^dx^dx" + if(p) [(1 + r 2 )<iA 2 + \ 2 dtf 2 \ . (2.6) 



Now the action for D6 brane becomes 



Sdg = -T 6 J d 7 ae- <t, ^/-det{g + 27ra'F), where T m = 1 . 
By using the well-known identity 



(2.7) 



det I ^ B j = det A ■ det{D - CA~ X B) (2.8) 



up to quadratic order in fields, we obtain 



£0 = "S^( 1 + SO l + \(F^F^ + 2F Xll F x n 

where F^ = d^Ay - d u A^. 



(2.9) 



3. Meson spectroscopy 

In this section we are to compute the vector meson masses by considering gauge field 
fluctuations on the D6 brane. Although we have two different scale of meson masses, 
light mesons and heavy quarkonia, the origin of the difference is that of the quark mass 
rather than the interaction. So we need to introduce only a single interaction scale Mkk- 
This scale is to be matched with the gauge theory scale Aqcd, hence it is encoded in the 
background geometry while the quark masses are encoded in the geometry of the probe 
branes. 
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3.1 Embedding 

We first find a D6 embedding geometry by solving the equation of motion for r(A). From 



the DBI action ( |2.9|) , we obtain the equation of motion for r(A) to be 
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With the following dimensionless variables 

A — > Uk\ t — )• Ukt, p^-UkP- 

we rewrite it as 
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(3.1) 



(3.2) 



(3.3) 



For large A, we can solve the equation of motion for r(A) to obtain the asymptotic solution 



as 



r(X) ~ + 



A 



(3.4) 



where Too, c are related to the quark mass and the chiral condensate, see j7| for details. 
3.2 Scalar and pseudo-scalar fluctuations 

We start with scalar and pseudo-scalar fluctuations. Though these were extensively studied 
in 0, we include them for completeness, not to improve the results in fiffl. The transverse 
fluctuation of the D6 brane is given by 



r(a^, A) = r„(A) + <$r(z"), 0(x'\ A) = 6<f>(x», A) 



(3.5) 



where r v (X) is the solution of the embedding equation. Inserting ( |3.5| ) into the induced 
metric (p.6j) and the DBI action (|2.9D, we obtain the induced metric 



ds 2 



— J i lllv dx^dx u + K[(l + r 2 )ciA 2 + A 2 dO|] + 2Kr v d a 5rdXdx a 



3/2 



+ J fT[d a <5rd 6 5r + rld a 5(/)d b 5(f)]dx a dx b , 



(3.6) 



where a and 6 run over to A and the DBI action, up to quadratic order, 
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Now we arrive at the linearized equation of motion for pseudo-scalar, 



= dx 



X 2 r 2 
V 1 Tr* 



r2A 2 yTT^ rr2 9M| 



4Mf 



(3.8) 
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and for scalar, 



= d x 



X 2 



;i + ^)3/ 2 v 1 + v' dx6r 



(3.9) 



X 2 U 2 9 M 2 r 
U 3r -\ 5r-d x ( ^ Tvfv 5 ( I + — ! | 



1 

4pf 



These equations are numerically solved to get meson masses with proper boundary condi- 
tions 0. 

3.3 Gauge field fluctuations 

Now we move on to the gauge field fluctuation. The relevant part of the Lagrangian density 
for the gauge field is given by 



^VA 2 ^/^+^^V /a f/ r -V / V^ F | 2 *TF*F vX 1 



1 + 

where is the embedding solution. We decompose the gauge fields in terms of the or- 
thonormal basis ip n , 4> n as 

A^,X) = J]bW(x%(A), A a (/,A) = J>("V)0 n (A) . (3.11) 

n n 

Then, the field strength takes the following form 

F, u {x^ A) = ^fW(/)i(A), FW(^) = - , 

n 

F^\(x^, A) = ^ tap (n) (aO0»(A) - c^ n (A)B<f A . 

With the decomposition, the quadratic part of the Lagrangian for the field reads 

R\ ^ _c m i,,„ , , 2 



(3.12) 



C 



B 



9p 
4 



E 



K(p)(l+f2 



(3.13) 



where = (l + A 2 -y/l + (p 7 ) 3 ^- To recover the canonical kinetic term of the 



gauge field in 4D, we impose the normalization condition for the wave function ^(A) as 



(2-Ka ) T 6 / dXy/go [ — ) ifj m ipn = fir, 



R 



3/2 



(3.14) 



where T§ = —T^V^/gs and V2 = J d^yjgmgn- We will impose a similar condition for cj>(\). 
The wave function tl>(X) satisfies the following mode equation derived from the quadratic 
action 



K(p)(l + r 2 ) 



9n 1 J] 1 J 



(3.15) 
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where = — q 2 . Then, we obtain 

(2na') 2 n [ dX^ 



K(p)(l+r 2 ) 

where m n is the eigenvalue. From ( |3. 14 ) and ( 3.16 ), we have now 

S m = -N I d 4 xJ2 (-F^F^ + \m\B^B^ 
J n=l ^ 

We rescale the coordinate by Uk to obtain 



go y m; 



U 3 / 2 4 M 



(3.16) 



(3.17) 



(3.18) 



KK 



To solve this equation, we impose two boundary conditions: at the IR, either VVi(O) or 
■0n(O) = 0) an d at the UV, ip n ~ X a with a < 1/2 from the normalizability condition in 

dm. 

Now we consider 4> n . The normalization of n is similar to ( 3.16| ) 

1 



(2na') 2 n I dXy^Q- 



] K{p){l + r 2 ) 



i0n — &n 



(3.19) 



As explained below and also in [|T4[, once we set the field as 4> n = ijj n /m n for n>l, it can be 
gauged away as a part of field. However, zero mode is exceptional, which is orthogonal 
to the other modes 



m r . 



T& / dXJTfo 



1 



K(p)(l + r 2 ) 



rtn = (forn>l) 



If we take O = CK (p)(l + r 2 )/^, 

(00, 4> n ) = dXip n = V'n(oo) - ip n (0) = (for n > 1). 

J 

Then, the constant C is given by 

1 = (0o,0o) -> C= ((2W) 2 f 6 | ^ A>) ^ + ^ )' 1/2 - 
The field strength is written as 

F^Or", A) = ^°(*)0 O (A) + ]T ( m -%^ n \x) - B«] ^ n (A) 

By gauge transformation, 2?^ absorbs d^ip^, 

and therefore the action ( |3.13| ) becomes 

S D6 = f d 4 x \d^d^ + f; (-F^F^ + i^BWsW" 

L n=l ^ 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 
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Note that for the heavy quarkonium system, 

-2/3\ " X /2 



o=( ( w)^ J f«^3(i + ^ 

/ f 00 1 \ _1 / 2 / 1 

= ( (2 7 ra , ) 2 T 6 jf dA^J = [ (2vra , ) 2 T 6 - 



-1/2 

= 0, (3.26) 



and so (fo =0 due to ( |3.22| ). To impose the UV boundary condition for tp n more precisely, 
we consider the mode equation ( |3.18| ) at large A, 

dx(\ 2 d x ip n ) = -^ n . (3.27) 

With ip n ~ A°, we obtain o;(ct — 3) — 0. Since the normalizability condition dictates 
a < 1/2, we should choose a = 0. 

3.4 Numerical results 

We solve the mode equations for scalar, pseudo-scalar, and gauge field fluctuations numeri- 
cally. We first compute the mass of light mesons to fix the model parameter and M^k- 
Since D4/D6 model has no non-Abelian chiral symmetry but for U(l)x, the pseudo-scalar 
meson in this model corresponds to rj' in QCD [jTj. In QCD, however, U(1)a symmetry is 
explicitly broken by the axial anomaly, and the observed mass of r/, my = 958 MeV, is 
much larger compared to the pion or kaon mass. Note that some portion of the r( mass 
comes from the anomaly effect which scales as Nf/N c . Since we are working in the large 
N c limit, we may use the mass of rj' with the anomaly contribution turned-off. So we use 
non-anomalous rf mass to obtain a rough number for the model input. To this end, we use 
the mass relation for Goldstone boson obtained in chiral perturbation theory at large iV c 



MY- 



2m„£ 2£(2m„ + m s ) 6r 
ml = — §-, m 2 ,= { g + (3.28) 

J IT ^Jn J TV 

where m u = rrid = m q . The term with r is from the axial anomaly. Now we take r = to 
estimate the rj' mass from non-anomalous contribution. With m q = 7 MeV, m s = 150 MeV, 
f w = 93 MeV, and S = (230 MeV) 3 , we obtain ~ 140 MeV and m' v ~ 390 MeV. Note 
that the mass of qq bound state such as p meson mass stays almost constant as we increase 
N c : for instance the light meson mass at large iV c is extensively studied in a unitarized 
chiral perturbation theory |l7]] and also in lattice QCD [jig] . 

We use /3-meson mass and rj mass in large N c limit as inputs to fix r^, = 0.191 and 
Mkk = 1-039. Our fitting results are summarized in Table 2. In Fig. ffl, we plot the masses 
of scalar, pseudo-scalar and vector as a function of r m . As observed in |7j , for large 
the meson mass becomes degenerate and increases monotonically with r^. This is simply 
because the equation of motion for scalar, pseudo-scalar and vector for heavy quark system, 
roo ^> Uk, are degenerated. 

L 3 Af 2 1 
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Figure 1: Meson masses for scalar (blue), pseudo-scalar (red), vector mode (dashed). The pseudo- 
scalar meson mass vanishes when current quark mass goes zero, see red line. For the high enough 
quark mass, all the meson masses are degenerated because of supersymmetry restoration 0. 



where can be real scalar 5r, pseudo-scalar (f), or vector tp, and y is a rescaled coordinate 
y=A/r oc . In the equation ( 3.2S| ) , the only scale is 

L 3 m„M K 



M 2 "" q ' " lx . (3.30) 

This means that for the heavy quark system, the fluctuating field has mass being propor- 
tional to \J mq \ t lK ■ As discussed in Jjj, the reason for the degeneracy is supersymmetry 
restoration. They are all in the same supermultiplet and for small m q limit supersymme- 
try is broken so their masses split. However, for the large quark mass or large separation 
between D4 and D6, the embedding is nearly flat, and so D6 brane restores supersymmetry. 



Mode 


Input (MeV) 


M/M KK 


M (MeV) 


Ps (rf) 


390 


0.375 


390 


Rs (<r) 




0.918 


954 


V (p) 


770 


0.741 


770 



Table 2: Light meson masses to fix free parameters, Mkk an d Too, in the model. 



Now we move on to the heavy quarkonium. To describe the light meson and heavy 
quarkonium systems on the same footing, we use the value of Mkk determined in the light 
quark sector to calculate the mass of heavy quarkonium. We choose J /ip mass as an input 
to fix r^. In Table 3, we list our results on the charmonium with = 11.92 , and in 
Table 4 we show bottomonium masses = 108.15. Note that r m is related to the quark 
mass parameter of D4/D6 model as 1 

m « - w - = < 3 - 31 > 



1 As well-known, the quark mass in D4/D6 model could be different from that in QCD by a constant 
factor. To obtain the constant we need to compare the scalar two point function obtained in D4/D6 model 
with that in the operator product expansion of QCD, for example see Jlj] ]. 
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Mode 


Particle data 


M/M KK M(GeV) Error 


Ps (n c ) 
Re ( x ) 

v (J/V) 


2.980 GeV 
3.414 GeV 
3.096 GeV 


2.978 3.095 3.72 % 

2.979 3.096 9.31 % 
2.978 3.095 (*) 



Table 3: Charmonium mass. 



So far we have tried to study light mesons and heavy quarkonia on the same footing, i.e. 



Mode 


Particle data 


M/M KK M(GeV) Error 


Ps (Vb) 

Rs (x) 

V (T) 


9.389 GeV 
9.859 GeV 
9.460 GeV 


9.105 9.46 0.75 % 
9.105 9.46 4.05 % 
9.105 9.46 (*) 



Table 4: Bottomonium mass. 



using the same value of Mkk for both cases. As long as the lowest lying KK mode masses 
are concerned, this unified approach seems working apart from the axial-vector meson. 
Note that this defect of degenerate vector and axial- vector mesons is the feature of D4/D6 
model itself not that of unified approach. However, the way this attempted unified approach 
works seems almost guaranteed in the following sense. The light meson masses are mostly 
inputs to fix the and Mkk, while for heavy quarkonia r 1 ^ will be fixed to reproduce a 
heavy quarkonium mass. As in Fig. (|l]), the masses of heavy quarkonium (r^ = 10.72) are 
almost degenerate regardless of their quantum numbers, while experimentally those masses 
are not very different from each other. Therefore, our results on the heavy quarkonium 
will be within experiments by roughly 10% deviations. 



4. Heavy quarkonium in deconfined phase 

We consider the meson spectroscopy in deconfining phase. By double Wick rotation, we 
get the black D4 background which is dual to the deconfining phase, 

ds 2 = (^) 3/2 (-f(U)dt 2 + dx 2 + dr 2 ) + ^ + U 2 dn 2 

^ = 9s(j) 3/ \ F i = dC 3 = ^e 4 , f(U) = l-^. (4.1) 

Here the temperature of the black D4 brane is identified with the temperature of a dual 
gauge theory 

6t ° =1 r Sr = lt K - (4 - 2) 

In terms of these periods, the metric parameters are given by 
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As discussed in M, there exists first order phase transition between confining and black 
D4 brane backgrounds; this phase transition is a deconfinement phase transition in field 
theory side. The critical temperature of the transition is fixed by the condition 5r = 5tE 



- dec 



Mk 
2vr 



3 

4-7T 



Uk 
L3 • 



Note that the temperature of the gauge theory is given as 

M . / 9 m q M KK 



_3_ 

4tt 




where M 



4vr X t 
4.1 Embedding 

The D6 brane embedding in the black D4 background is described by r(A) 

'[A 3/2 

The DBI action for r(A) is then 

1 



(4.4) 



(4.5) 



ds 2 = ^ 



f(U)dt 2 + dx 2 



+ K(p) 



;i + r 2 )d\ 2 + x 2 dn 2 



(4.6) 



Sd6 



d ae 



= -T m J d 7 aVh\ 2 ^f(U)(l + ^) (l + ^ 

From this, we obtain the equation of motion for the embedding function r(A) 

X 2 r 
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16^6 



vT+7 5 



0. 



(4.7) 
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The solution of the embedding equation was extensively studied in |JJ. By comparing the 
energy density of various embedding solutions, the authors of 0] found that there is a 
first order phase transition at ~ Uk- This transition, from Minkowski to black hole 
embedding, occurs when the D6 brane touches the black hole horizon at sufficiently large 
temperature. The transition temperature for charm quark is given by Tf unc i ~ 1.0202M C , 
where M 2 = j^M KK = (0.571GeF) 2 . With the parameters fixed in confining phase, we 
obtain T^ ec ~ 165 MeV and Tf un d ~ 582 MeV. We summarize the different embedding 
solutions below. 



T 


Embedding 


Energy spectrum 


T < Tdec 

T dec < T < T f^nd 

Tf U nd < T 


Confining emb. 
Minkowski emb. 
Blackhole emb. 


discrete 
discrete 
continuous 



Table 5: Classification of embcddings 



We note, however, that there is no priori reason for T^ ec < Tf un d, which is satisfied 
only when 

II OS 

(4.9) 



0.98 
~9tT 



m f. s . > z^\ t M KK ~ 36 A t MeV. 
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As a consequence, there is no Minkowski embedding for the light quark system, meaning 
that light mesons will be melted away immediately when the temperature exceeds the 
deconfmement temperature T > Td ec - For the heavy quark system we expect that there 
are discrete meson spectra even in the deconfining phase since is big enough to satisfy 
the condition. 



4.2 Fluctuating fields as mesons 

In this section we study temperature dependent meson masses. The temperature depen- 
dence of scalar and pseudo-scalar heavy quarkonia are already considered in ||] , and now we 
extend the work by including vector mesons. In the high-temperature deconfining phase, 
there is only one relevant scale, T/M, which is nothing but the separation between D4 and 
D6 branes. Therefore varying temperature is equivalent to varying the D4-D6 separation 
with a fixed quark mass. 



4.2.1 Scalar and pseudo-scalar fluctuations 



Scalar and pseudo-scalar mesons are described by the D6 brane fluctuating on the trans- 
verse direction, r(A) = r v (X) + 5r(t,X), 4>(t,\). The quadratic part of Lagrangian after 
rescaling by Ut is 



s/h 



2^1+T 2 
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16p%J 1 + r 2 



(d x 5rf + 37^ 5r2) + / 3(1 + r 2 ) (A 2 - 7r 2 ) 
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1/3 4(l + r 2 K 2 



1 



Ap%) l-Apl 

The equation of motion for the bulk scalar field (<5r) and for pseudo-scalar (0) 
1 \ X 2 d x 5r 



= d 
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1 \V3 



l&plJ (l + r 2 )3/2 
4A 2 



3A 2 (A 2 -7r 2 )yi+7 2 3 



,io 







A 2 r„r„ 



8 \Pv y/l + i* 



Sr 



v / T+^(i-4p§) 



1 + 
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w 2 



1 



i-^) *| ^ (4.11) 



where w 2 = {f^w, k 2 = {f^k 2 , and the equation of motion for the pseudo-scalar is 
1 \ A 2 r 2 



= dx 



1 



4p? 



kl u. 



(4.12) 

Since there is no Lorentz symmetry due to the temperature, we define the thermal meson 
mass as d 2 cf) = m\(j), i.e., at vanishing spatial momentum. The fields 4>, 5r have two linearly 
independent solutions at A — > oo as A , A -1 . Here we choose the asymptotic solution as 
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Figure 2: The temperature dependent meson masses. Real scalar (left) and pseudo-scalar (right) 
masses in GeV unit. In confined phase the temperature dependence is trivially constant which is 
smooth at the deconfinement transition temperature, Tdec ~ 160 MeV. Note that the disjoint at Tdec 
denotes that we use different backgrounds below and above Tdec- Above the melting temperature, 
Tfund ~ 571 MeV, there is no discrete spectrum for the meson. 




Figure 3: Left : Dispersion relation for the real scalar with various temperature, T = (1.15, 2.6, 
3, 3.24) Tdec from top to bottom. Middle: J as a function of k. Right: The group velocity 
v g = ^|fc^oo as a function of temperature. 



sion relation of heavy quarkonia, which is important to understand the dissociation and 
the screening mass of heavy quarkonium. The result in Fig. || shows that for sufficiently 
low temperature but larger than Td ec , the group velocity v g is approaching one, and it is 
decreasing with increasing temperature. Note that in D3/D7 model the dispersion relations 



of mesons was studied [20|. 



4.2.2 Gauge field fluctuations 

Now we move on to the gauge field to study vector mass in deconfined phase. The relevant 
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Lagrangian is 



Co = ^VTT^ (i - (f^ + u%f*>) 



"Pg^F^Fap + 2g xx g^F Xfl F Xa 



(4.14) 



Transverse vector 

The vector fluctuations at finite temperature are classified into transverse and longitudinal 
modes due to lack of Lorentz invariance. We first consider the transverse part. Assuming 
that the boundary gauge field propagates only in x direction, we obtain the Lagrangian for 
transverse vector fluctuations as 

r 



L G = 2V50 



1 



(d t Ayf + {d x A y f + 



K{ 



f(U) 
tran 



i 



(4.15) 



The equation of motion for the transverse vector is 

3/2 



dl + dl)A y 



f 



0, 



(4.16) 



where ^/go = A 2 Vl + r 2 



u' T 

4pJ 



3 \ 2 



(f) 3//2, ^ e ca l cu l a t e the dispersion relation and 



v,(J/¥) 
l.Or 




Figure 4: Dispersion relation for the transverse vector with T = 1.04, 2.91, 3.37 Tdec from bottom. 

group velocity to obtain Fig. [|. Note that w ~ k 2 for small k and w ~ k in large k region. 
Longitudinal vector 

We decompose the gauge field into the ortho-normal basis 



to obtain 



C G 





A; 






4 


[(* 


Vdo 




4 


[(* 



(4.17) 



3/2 



7 / 



/ 



K(l + r 2 

ib~.il. 

+ 2 



if(l+r 2 ) 1 



(4.18) 
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Figure 5: The temperature dependent vector meson masses in GeV unit. Note that the disjoint 
at Tdec denotes that we use different backgrounds below and above Td ec - 



Again, we define the thermal meson mass at zero momentum, <9j = or hi = 0. Then the 
wave-function of the gauge field satisfies 



d\ 2Y 



3/2 

u) kVT+^ 



(4.19) 



with the normalization condition 



r \2 



(4.20) 



From these, we obtain the equation of motion for the vector mesons 




3/2 




where My = j^rriy is dimensionless. The meson mass is given by 



V i>m = 0. 



(4.21) 



my = \ —^My- 

V L 6 



(4.22) 



We solve the mode equation numerically and show the result in Fig. B We impose 
tp m ~ 1/ A at large A due to the normalizability. 

As shown in Fig. | and Fig. g the mass changes smoothly from confining to deconfining 
phase. This may be different from previous studies based on a bottom-up AdS/QCD 
model H, QCD sum rule approach [21], and holographic heavy potential p2| . Note, 
however, that the maximum shift of charmonium mass at T = 1.05T C is about 200 MeV pifl, 



which is roughly 6.7% of the mass of J/t/j and may be too small to be explained by our 
approach based on large iV c approximation. 
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4.3 Remarks 

Now, we discuss the relevance of the D4/D6 model for light mesons and heavy quarkonia. 
One of our basic goals is to study both light mesons and heavy quarkonia on the same 
footing based on a single D-brane model with a common energy scale Mkk- Another 
one is to study heavy quarkonium properties in the high-temperature deconfining phase. 
Apart from non-Abelian chiral symmetry, the D4/D6 model might be a good candidate 
for these purposes since it has both confining and deconfining phases and includes quark 
mass through embeddings naturally. Especially for a heavy quark system, the D4/D6 
model shows some similarity with heavy quarkonium physics. First, the non-Abelian chiral 
symmetry is not an issue for heavy quarks. Second, the D4/D6 model has one more 
transition at T = Tf un d other than the deconfinement transition, which may be associated 
with the dissociation of heavy quarkonia at high temperature. A positive clue for this is 
that Tf un d is proportional to the square root of a quark, meaning that the charmonium 
melts relatively at low temperature compared to the bottomonium. For instance, the 
dissociation temperature for a bottomonium T is ~ 2.7T C , while that for a charmonium 
J/V is ~ 1.3T C |§. 

We finish this section with a summary of the discussion in l| on the usefulness of 
Dq/Dp systems in studying meson bound states. A Dq/Dp system may be good for 
ss bound states at high temperature since the mesons in the Dq/Dp system are deeply 
bounded. Even above Tf un d, there still exists some broad peak in the spectral function of 
two-point correlators of mesons. Therefore Tf un d is slightly smaller than the dissociation 
temperature of the heavy quarkonium. The D4/D6 model may describe some exotic gauge 
theories. However, there exist certain properties of heavy quarkonia in the quark-gluon 
plasma that could be understood in the D4/D6 model. 



5. Summary 

We re-analyzed the D4/D6 model to study holographic light vector mesons and the prop- 
erties of heavy quarkonium in confining and deconfining phases. To treat the light mesons 
and heavy quarkonium on the same footing, we used the same compactification scale Mkk 
in both systems. In confined phase, we observed that the meson spectroscopy of the light 
meson and heavy quarkonium could be described with a single Mkk- We found that 
like scalar and pseudo-scalar mesons the vector meson mass is linearly proportional to 
the square root of the quark mass, when the quark mass is large. With a Mkk fixed by 
light meson masses, we calculated the mass of heavy quarkonium in confining and decon- 
fining phases. Unlike previous studies based on a bottom-up AdS/QCD model ||, QCD 



sum rule approach [21|, and holographic heavy potential |22], we found that the mass of 



heavy quarkonium changes smoothly from confining to deconfining phase. However, since 



the maximum shift of charmonium mass at T = 1.05T C is about 200 MeV [21], which is 
roughly 6.7% of the mass of J/t/j, it is not conclusive if our results are really different from 
a previous study. We also obtained the in-medium dispersion relation for heavy quarko- 
nium, which is important to understand the dissociation and the screening mass of heavy 
quarkonium. 
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Certainly there are many things to be improved in our study to be more close to 
QCD. We list some of them here. Surely chiral symmetry should be the first one. As 
well known, non-Abelian chiral symmetry is essential to understand light mesons, and it 
is also important for heavy-light system due to the light quark. Second thing is how to 
include heavy-light meson in this picture with correct chiral symmetry and heavy quark 
spin symmetry. 
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